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Abstract 

Results for the first test of the "crude" QCD continuum model, commonly 
used in QCD Sum Rule analyses, are presented for baryon correlation func- 
tions. The QCD continuum model is found to effectively account for excited 
state contributions to the short-time regime of two-point correlation functions 
and allows the isolation of ground state properties. Confusion in the literature 
surrounding the physics represented in point-to-point correlation functions is 
also addressed. These results justify the use of the "crude" QCD continuum 

model and lend credence to the results of rigorous QCD Sum Rule analyses. 
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I. INTRODUCTION 



Point-to-point correlation functions have been the subject of intense study since the 
early days of lattice field theory. In the quest for an ah initio determination of the low- 
lying hadron spectrum, lattice QCD investigations have focused on the large Euclidean-time 
tails of three-momentum-projected two-point correlation functions. As such, the short time 
regime, where excited state contributions are significant, has simply been discarded. 

The lattice approach to QCD allows the determination of correlation functions deep in the 
nonperturbative regime. In contrast, the QCD Sum Rule (QCD-SR) method [|I[] is restricted 
to the near perturbative regime of the truncated Operator Product Expansion (OPE). In 
this regime, one cannot ignore the contributions of excited states in a phenomenological 
description of QCD correlation functions. To account for these contributions, the so-called 
"crude" QCD continuum model is introduced |]]. This model exploits the leading terms of 
the OPE, and introduces a sharp threshold marking the onset of the QCD continuum. The 
contributions of this model relative to the ground state, whose properties one is really trying 
to determine, are not small. They are typically 10 to 50% 0. The validity of this model is 
relied upon to cleanly remove the excited state contaminations. 

This investigation examines the physics in the near perturbative regime of point-to-point 
correlation functions where QCD-SR analyses are performed. The QCD continuum model 
is constructed for three- momentum-projected Euclidean-time two-point functions following 
the techniques established in the QCD-SR approach. This model will then be used as a 
probe of the physics represented in lattice QCD correlation functions, and as a test of QCD- 
SR techniques. This investigation explores some of the ideas briefly summarized in Ref. 
in greater depth and detail. 

Some attention has recently been given to the behavior of lattice point-to-point corre- 
lation functions in the near perturbative regime [§]||. There the focus was on correlation 
functions in coordinate space. It was concluded that the QCD-SR-inspired continuum model 
was sufficient to describe the lattice correlation functions over the calculated range. Using 
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the entire lattice correlation function (including the deep nonperturbative regime), ground 
state masses were extracted and were found to agree with conventional lattice analyses. 

Here the emphasis is on determining ground state properties by examining only the 
first few points of the lattice correlation function in the near perturbative regime. This is 
similar in spirit to QCD-SR analyses. In this manner, the validity of the QCD continuum 
model is rigorously tested. By extending the analysis interval of the correlation function 
deeper into the nonperturbative regime, the evolution of fitted ground state properties may 
be monitored. A sensitivity to the analysis interval in the extracted parameters would 
indicate a failure of the QCD-SR-inspired continuum model. Comparisons are made with 
conventional lattice results where ground state properties are simply extracted from the 
tails of the correlation functions. The importance of higher order terms of the OPE in 
the formulation of the QCD continuum model is also examined. These terms were not 
investigated in Ref. |H|. 

One would also like to evaluate whether these techniques are useful in analyzing other 
lattice QCD correlation functions. Such techniques may be of practical use for analyzing 
two-point correlation functions which become noisy prior to a clear ground state domination, 
such as heavy-light meson correlators. In a subsequent paper 0, these techniques will be 
exploited to investigate nucleon properties obtained from unconventional interpolating fields. 

The outline of this paper is as follows. Section |I| introduces the definition of the two- 
point function examined in this investigation and summarizes the lattice techniques and 
parameters. In Section |T|, the QCD continuum model is derived for three-momentum- 
projected two-point correlation functions in Euclidean space. Issues associated with relating 
the lattice and continuum (lattice spacing a — > 0) formulations of the model are discussed 
here. Fits of the nucleon correlator are presented in Section [TV]. The validity of the QCD 
continuum model as implemented in the QCD Sum Rule approach is evaluated in Section 
|V[ The fit parameters are compared with other approaches in Section [VI]. In section |V1I| , 



the physics represented in two-point correlation functions is discussed, where erroneous 



conclusions in the literature are addressed. Finally, Section [VIII] summarizes the findings of 



this investigation. 



II. LATTICE CORRELATION FUNCTIONS 

A. Interpolating Fields 

In lattice calculations, the commonly used interpolating field for the proton has the form 

X i(x) = e abc (u Ta (x)C l5 d b (x)) u c (x) . (2.1) 

Here, we follow the notation of Sakurai |7]]. The Dirac gamma matrices are Hermitian and 
satisfy {7^,7^} = 2^, with a^ v = ~ [7^,7^]. C = 7472 is the charge conjugation matrix, 
a, b, c are color indices, u(x) is a w-quark field, and the superscript T denotes transpose. 
Dirac indices have been suppressed. 

In the QCD Sum Rule approach, it is common to find linear combinations of this inter- 
polating field and 

X *(x) = e abc (u Ta (x)Cd b (x)) 75 « c (x) , (2.2) 

which vanishes in the nonrelativistic limit. With the use of Fierz relations, the combination 
of the above two interpolating fields with a relative minus sign may be written 

Xsr(x) = e abc (u Ta (x)C lllU b (x)) l5 rd c (x) , 

= 2 (X2 - Xi) , (2-3) 

giving the proton interpolating field often found in QCD Sum Rule calculations |2|]8|.[| • 

Here the task is to evaluate the validity of the QCD-SR-inspired continuum model. This 
issue is best addressed with the use of the interpolating field Xi °f ([2-1D- For comparison 
with other calculations, results for xsr are a ls° reported. 

The A resonance is also considered for comparison with the nucleon results. For sim- 
plicity, the A ++ interpolating field is used, 

x £(x) = e abc (u Ta (x)C^u b (x)) u c {x) . (2.4) 
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B. Two-Point Function 



Two-point correlation functions are used to determine hadron masses. Consider the 
following two-point function for the nucleon, 



G 2 (t,p) = £ e -**tr[r 4 ( | T{xi(x)xi(0)} | ) 



(2.5) 



Here, T 4 = (l + 7 4 )/4 projects positive parity states for p — 0, and tr indicates the trace over 
Dirac indices. Correlation functions at the quark level are obtained through the standard 
procedure of contracting out time-ordered pairs of quark field operators, 



G 2 (f,# = 5>-**tr 



r 4 e aoc e 



abca'b'c' I nao' 



Sf(x, 0) tr Sf (x, 0) CSf T (x, 0)0-' 



+ST (x,0)CSf T (x,0)C- i S? (x,0) 



-1 ebb' , 



(2.6) 



where C = C75, and Sf{x, 0) = T {u a (x), u a '{0)}, etc. 



C. Lattice Techniques 



Here we briefly summarize the lattice techniques used to calculate Q2.6p . Additional 
details may be found in Ref. [ID|. Wilson's formulation is used for both the gauge and 
fermionic action. S"?7(2)-isospin symmetry is enforced by equating the Wilson hopping pa- 
rameters k u = Kd = k. Three values of k are selected and are denoted K\ = 0.152, k 2 = 0.154 
and k 3 = 0.156. To make contact with the physical world, the mass, continuum threshold 
and interpolating field coupling strengths calculated at the three values of k are linearly 
extrapolated to k ct = 0.159 8(2) where an extrapolation of the squared pion mass vanishes. 
Differences between linear extrapolations to = as opposed to the physical pion mass 
are small and are neglected in the following. 

Twenty-eight quenched gauge configurations were generated |10| by the Cabibbo- 
Marinari |TT] pseudo-heat-bath method on a 24 x 12 x 12 x 24 periodic lattice at (3 — 5.9. 
Configurations were selected after 5000 thermalization sweeps from a cold start, and every 



1000 sweeps thereafter [12 



Dirichlet boundary conditions are used for fermions in the time direction. Time slices 
are labeled from 1 to 24, with the 5-function source at t — 4. To minimize noise in the 
Green functions, the parity symmetry of the correlation functions, and the equal weighting 
of {[/} and {U*} gauge configurations in the lattice action are exploited. The nucleon mass 
determined from xi of ( |2.1| ) is used to set the lattice spacing. This estimate lies between 
other estimates based on the string tension or the p-meson mass. The lattice spacing is 
determined to be a = 0.132(4) fm, and a" 1 = 1.49(5) GeV. 

Statistical uncertainties in the lattice correlation functions are estimated by a single elim- 
ination jackknife ||13|| . A covariance matrix fit of the pole plus QCD continuum model over 



a typical range of 15 time slices is likely to be unreliable for 28 gauge configurations |14 



Instead we use standard statistical error analysis in which correlations among the fit param- 
eters are accounted for. The Gauss-Newton method is used to minimize \ 2 - Uncertainties 
are taken from the standard error ellipse |L5| at x 2 — Xmm + 1- 



D. Mean-Field Improvement 



In this analysis, the nucleon coupling strength (or residue of the nucleon pole) is deter- 
mined in absolute terms, without resorting to a ratio of the QCD continuum contributions 
as done in 



Our approach requires the use of mean-field improvement |lq|L7 

/2k 



4> 



Continuum 



7 3 / 2 



■0 



Lattice 



(2.7) 



where 



2k 



3k \V3 



4:K r 



(2.8) 



to account for otherwise large renormalization factors. While ( |2.8[ ) is derived for heavy 
quarks |IBfl , its use in light quark studies has already been suggested in fTTH - The k depen- 
dence of these two wave function normalizations is very different and use of the mean-field 
improved normalization is crucial to recovering the correct mass independence of the Wilson 
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coefficient of the identity operator. This point will be further illustrated in the following 
section. 

The implementation of Wilson fermions on the lattice induces mixing between the com- 



posite nucleon interpolating fields of ( [2.1[ ) and (2.2), reflecting the breaking of chiral sym- 
metry [Tj|. This mixing is discussed in more detail in Ref. The important point is that 
the mixing of the interpolating fields xi an d X2 is negligible. As a result, it is possible to 
identify the properties of these interpolating fields determined on the lattice with those of 
their continuum (a — > 0) counterparts to a good approximation. The principle renormal- 
ization constant C\ has been determined in the mean-field approach |T?| and is used in the 

following. The renormalization at the scale of 1/a is 

7 , o K N 3/2 

„ Continuum Xn ..Lattice M l fo n\ 

X = ^ X l 1_ 4^J ' (2 ' 9) 

and Z XN = (1 - 0.73 a v ) ~ 0.80 at (3 = 5.9. 



III. THE EUCLIDEAN SPACE FORMULATION OF THE QCD CONTINUUM 

MODEL 

A. Spectral Representation 

At the phenomenological level, the two-point function is calculated by inserting a com- 
plete set of states N l between the interpolators of ( |2.5| ) and defining 

(0\Xi(0)\N\ P ,s)=X\u( P ,s). (3.1) 

Here, the coupling strength, X\, measures the ability of the interpolating field xi to annihilate 
the i'th nucleon excitation. For p = and Euclidean time t — > 00, the ground state dominates 
and 

G 2 (t) -> Xle- AlNt . (3.2) 

We note here that the exponential suppression of excited states is somewhat similar to that 
encountered in Borel improved QCD-SR analyses. There the square of the mass appears 
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in the exponential as exp(— M^r 2 ) where r is the inverse Borel parameter. The spectral 
representation is defined by 

/■OO 

G 2 (t) = / p(s)e- st d Sj (3.3) 
Jo 

and the spectral density is 

p(s) = \j5(s-M N ) + as), (3.4) 

where £(s) provides the excited state contributions. 

While it may be tempting to fit the correlation function in the shorter time regime by 
including additional poles in the spectral density, 

p(s) = X 2 N 5(a - M N ) + \ 2 N * 5(a — M N * ) + • • • , (3.5) 

such an approach fails in a number of ways. If one simply includes a couple of poles for the 
spectral density, one finds that the mass of the excitation is completely determined by the 
leading time slice that is considered in the fit. This is illustrated in Figure [TJ where a fit of 
time slices 7 through 20 is shown. If time slice 6 is included, the mass of the second pole 
increases to accommodate the earlier time slice with smaller uncertainties. 

One also finds that the coupling strength A at. exceeds Xn by an order of magnitude. 
Our physical intuition suggests the opposite. The coupling strength should decrease for 
excited states as these states are expected to have broader wave functions. The probability 
of finding three quarks at the same space time point is therefore smaller. In addition, 
such an approach is a poor representation of the physics believed to be there as it neglects 
the natural widths of the resonances and the QCD continuum of multiple hadron states. 
The correlation function is probably best described by many states of diminishing coupling 
strengths and increasing widths. It may be that the QCD-SR-inspired continuum model is 
an efficient way of characterizing this physics. 
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B. Operator Product Expansion (OPE) 



The form of the spectral density used in the QCD continuum model is determined by 
the leading terms of the OPE surviving in the limit t — > 0. The calculation of these terms 
proceeds in the usual fashion of contracting out time-ordered pairs of quark field operators. 
For the nucleon interpolating field of one has 

G 2 (t,p = 0) = £ tr [ T 4 e abc e a ' b ' c ' (sf tr (S b u b ' l5 C Sj^CTS) + Sf l5 C Sj bb ' C' 1 ^ S?) 

(3.6) 



In Euclidean space and coordinate gauge the quark propagator has the expansion 



<W = J_ roa' , 1 

q ~ 2vr 2 x 4 (2vr) 2 x 2 



1 £aa! 



2 2 3 



(: qq :) 5 aa ' + 



and (j2(t) has the following OPE 

G 2 (t) ~ 



3 ■ 5 2 / 1 28 m q a 14 m 2 a 2 
2V \¥ + 2E~t^ + 2E~t r ~ 
567T 2 (:qq:) a 3 



75 t 3 



+ 



The quark mass used in the OPE is obtained from the standard relation 



m a a = 

q 2 



1 /l 1 



(3.7) 



(3.8) 



(3.9) 



and is renormalized to one loop Hl9|,^0 



C. QCD Continuum Model Contributions 

The spectral density used in the QCD continuum model is defined by equating ( |3.3j ) and 
(PD- Thus, p(s) is given by the following Laplace transform pairs 



1 1 n n = 1, 2, 3, . . . 

— — > -. — - — rr s n , . (3.10) 

Ret>0 



The QCD continuum model is denned through the introduction of a threshold which marks 
the effective onset of excited states in the spectral density. This model is motivated by 
duality arguments as realized in the Q 2 dependence of the experimental cross section for 
e + e~ — ► hadrons. The QCD continuum model contribution is 

poo 6 1 c fc 

/ p(s) e~ st ds = e -»« E E TT^CnOn, (3.11) 

Js n=l k=0 K ' 1 

where C n and O n are the Wilson coefficient and normal ordered operator of the term t~ n in 
( p.8|) . The phenomenology of G2{t) is summarized by 

p(s)e- st ds, (3.12a) 

-0 

= A?e-^ + e|^e- ot ( (3.12b) 

\t 6 t 5 2 t 4 6 t 3 24 t 2 120 t J 

25 \ t 5 t 4 2 t 3 6 t 2 24 t J 

14 Tn 2 a, 2 flfo lfolfol 

25 \ t 4 t 3 2 t 2 6 t J 
567r 2 (gg)a 3 f 1 a 1 ag| \ 
75 \t 3 t 2 2 t J J ' 

Here we have introduced an additional parameter, ^, governing the strength of the QCD 
continuum model. Strictly speaking, £ = 1 in the continuum limit, a — > 0, but here is 
optimized with Ai, M^, and Sq to account for enhancement of the correlator in the short 
time regime due to lattice anisotropy. Its deviation from 1 is a reflection of the fact that 
we are matching the model formulated in the continuum limit to the lattice results at finite 
a. Ref. found the anisotropy to be large for x — xo < 6 for free quark correlators and 
remain large in their interacting simulation at (3 — 5.7. At (3 = 5.9 there is some hope that 
anisotropy issues will be less problematic for the Fourier transformed correlators presented 
here. However, at very short times the quarks are essentially free and the anisotropy must be 
accommodated. With this approach, the effects of lattice anisotropy are absorbed through 
a combination of a larger QCD continuum model strength (£ > 1) and marginally larger 
continuum threshold (sq). 
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D. Lattice Cutoffs 



Unlike coordinate space lattice analyses infrared lattice artifacts are not a significant 
problem for this approach. The Fourier transform weight exp(— ip ■ x) is correct for all 
propagator paths including those which wrap around the lattice spatial dimensions. 

The ultraviolet lattice cutoff may be modeled in a manner similar to that for the QCD 
continuum model. The upper limit of the integral over the spectral density of (3.12a) is cut 



off at A, and additional terms appear in ( |3.12b| ). As an example, consider the leading term 
of the OPE. In this case, 

I r°° 1 

* = h a 86 ^- (313) 

is modified to 



A 1 ■= „ , 1 At f 1 A 1 A 2 1 A 3 1 A 4 1 A 5 ) 

— s 5 e ~ st ds = -- e~ At { — + — + - — + - — + + } . 3.14 

o 5! t 6 \t G t 5 2 t 4 6 t 3 24 t 2 120 t J V ' 

Figure |2] illustrates the C$/t 6 term with the ultraviolet cutoff correction (dashed line) fit 
to the lattice data from t = 5 — > 7 by optimizing A and Cq. The corresponding curve for 
Ce/t 6 without the ultraviolet correction (solid line) is also indicated. The optimum value 
for A is 4.55(3) which is not too different from the momentum cutoff of tt. Of course, the 
two numbers need not agree since we are modeling the ultraviolet cutoff and matching a 
continuum formulation, (a — > 0), to the lattice. The correction term accounting for the 
ultraviolet lattice cutoff is negligible by t — t = 2. Rather than introducing an extra 
parameter in the fit, the correction is simply neglected, and all fits begin at t = 6 in the 
following. 

A similar modification of the correlator due to the lattice regularization was not ac- 
counted for in Ref. ||, when comparing hadron correlators calculated on the lattice with 
free quark propagators, and correlators using continuum, (a — > 0), propagators. Inclusion 
of these effects in Fig. 3 of Ref. would suppress the continuum curve at short distances 
such that it may continue to follow the diagonal elements of the lattice data more closely. 
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E. QCD Continuum Model Summary 



Figure |3] displays the total and individual contributions to the QCD continuum model 
term of (|3.12|) . For comparison, the total and individual contributions to the OPE of ( |3.8|) 
are illustrated in Figure [|. The higher order terms of the OPE become significant at large 
time separations. Therefore, a truncated OPE will have a corresponding upper limit in time 
separations within which reasonable convergence is maintained. In contrast, the exponential 
of the continuum threshold in the QCD continuum model acts to suppress the higher order 
terms. The relative contributions of the subsequent terms of the OPE to the QCD continuum 
model are well ordered throughout the fitting regime, and inclusion of the first few terms of 
the OPE is adequate. Indeed, the contributions from the m 2 and quark-condensate terms 
make negligible contributions relative to the identity and m q terms. These higher-order terms 
are not included in the following analysis. In summary, the phenomenological description of 
the correlation function is taken to be 

G 2 (t) = A?e- A ^ + e^e-^( (3.15) 

\ t e t 5 2 t 4 6 t 3 24 t 2 120 t J 

2877VZ fl_ £o 1 4 1 1 4X\ 

25 \t 5 t 4 2 t 3 6 t 2 24 t J J ' 

The shaded areas of Figure |] indicate the regions typically excluded form Borel- 
transformed QCD-SR analyses. The remaining region is far from the regime where the 
ground state pole dominates the correlation function. A rigorous uncertainty analysis is 



needed to ascertain the predictive power of QCD Sum Rules [21 . 



IV. NUCLEON CORRELATOR FITS 

Among the first of things to verify is the leading time dependence of G^if)- Here we 
consider the lattice data at n — 0.156 where m q corrections are smallest and one has the 
best chance to verify that the dependence is indeed proportional to 1/t 6 . Fitting time slices 
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t = 6 and t = 7 to the functional forms of t~ 7 , t~ G , and t~ 5 yields a x 2 /dof of 80, 4.8, and 
230 respectively, suggesting the predominant behavior is 1/t 6 as anticipated. 

Our purpose is to test whether the nucleon mass and coupling strength can be obtained 
accurately from a fit considering only the first few points of the correlation function. The 
lattice correlation functions |TD| are fit with ( p.l5| ) in a four-parameter search of Ax, Mat, 



Sq and £ in analysis intervals from t = 6 — > tf where tf ranges from 11 through 23. Figure 
H] illustrates these 13 fits of the lattice correlation function at the smallest value of k con- 
sidered. The statistical uncertainties are smallest for this quark mass and provides the best 
opportunity for revealing structure in the lattice correlation function that is not accounted 
for by the QCD continuum model. The typical x 2 /dof for these pole plus QCD continuum 
model fits is 0.6. 

The nucleon mass determined in each of the intervals is plotted as a function of tf in 
Figure |6|. Pole plus QCD continuum model fits are compared with simple pole fits. Figure 
|7| illustrates similar results for the coupling strength. It is interesting to see that the simple 
pole determination of Ai fails to form a plateau at large time separations. Similar results 
are seen for the larger values of k = 0.154 and 0.156. The fit parameters extracted for the 
region t = 6 — > 20 are summarized in Table | [^|. The plateau in M N and Ai for pole plus 
QCD continuum model fits indicates that the QCD continuum model effectively accounts 
for excited state contaminations in the correlation functions. The search parameters So and 
£ display a similar plateau as one might expect, since the regime in which these parameters 
are largely determined is common to all intervals. 

Table p] also indicates a value for £ obtained by fitting £ x the first two terms of the OPE 
to time slices t = 6 and 7. Higher order terms not included in ( |3.8|) are expected to be 
significant and therefore £ obtained in this manner is only a rough estimate. Its similarity 
to that obtained from the QCD continuum model verifies the deviation of £ from unity is a 
lattice artifact and not a failure of the QCD continuum model. It should also be mentioned 
that £ has been defined with respect to the Wilson coefficient of the identity operator only to 
leading order in perturbation theory. a s and leading log corrections are expected to reduce 
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£ by about 25%. 

A similar fit including the Wilson coefficient of the m q correction, C5, as a search parame- 
ter confirms the OPE value for the coefficient ratio C^/Cq. The fact that this estimate agrees 
with the OPE prediction confirms the role of £ as an overall continuum model strength. As 
such, £ is expected to be independent of the quark mass in the same manner that Cq, the 
coefficient of the identity operator, is. To demonstrate that this is in fact the case we present 
Figure || Here the quark mass dependence of £ is illustrated (solid curve) for the three values 
of k considered on the lattice as well as the value linearly extrapolated to k ct where the pion 
mass vanishes. £ is nearly independent of the quark mass. This behavior is also in accord 
with the fact that enhancement of the correlator due to lattice anisotropy is largest at very 
short time separations where the OPE is dominated by the identity operator. 

This result is not trivial, as it requires the use of the OPE value for C5 and the use of 
mean-field-improved operators. If the naive wave function normalization, v^2k, is used, £ 
varies as indicated by the dashed line plotted in Figure §. 

In summary, the pole plus QCD continuum model allows the extraction of nucleon ground 
state properties, even for an interval as small as t — 6 — > 11. The techniques examined here 
may be useful in analyzing correlation functions that suffer a loss of signal prior to clear 
ground state domination. These techniques are used in a subsequent paper investigating 
nucleon properties obtained from unconventional interpolating fields 0. 

V. QCD SUM RULE TEST 

In the sum rule approach £ is fixed to unity by the OPE. To test the QCD-SR method 
as closely as possible to its actual implementation, we will use the value for £ obtained in 
the previous fit from t = 6 —>■ 20. This value is similar to that obtained by fitting ( |3.8| ) to 
the first few time slices of the lattice data. 

Figure |9| illustrates the correlation function fits for the lightest quark mass at n = 0.156 
for the three-parameter search of Ai, and Sq, in analysis intervals from t = 6 — ► tf where 
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tf ranges from 9 through 23. While ideally we would like to start at tf = 7, a minimum of 4 
points is required for the fit with uncertainty estimates for the fit parameters. The typical 
X 2 /dof for these pole plus QCD continuum fits is 0.53. For the heaviest quark mass where 
the statistical uncertainties are smallest the typical x 2 /dof is 0.6. 

Using a pole plus QCD continuum model, one can extract ground state nucleon properties 
using only the first few points of the correlation function. The stability of the nucleon mass 



to the analysis interval is displayed in Figure [10]. The nucleon coupling strength is illustrated 



in Figure ITT] where it is plotted as a function of the analysis interval. All three fit parameters 
are stable as functions of tf. Moreover, the nucleon mass and coupling strength agree with 
the values extracted from the tail of the correlation function with a simple pole fit from 
tf — 7 — > tf when tf ~ 20 J23|. This indicates the QCD continuum model inspired by QCD 



Sum Rules is successful in quantitatively accounting for excited states in point-to-point 
correlation functions. However, interplay between the fit parameters gives rise to rather 
large uncertainties at tf — 9. 

The smallest analysis interval considered here is still somewhat generous relative to the 
interval considered in QCD-SR analyses. Borel improved Sum Rules including terms to 
dimension 8 in the OPE are typically analyzed in the regime r = 1/Meorei — 0.15 — > 0.35 fm 
corresponding to t ~ 5 — > 7 in the previous figures. However, Borel improved sum rum rules 
suppress excited states more effectively than the Euclidean formulation presented here, as 
the phenomenological side of Borel improved sum rules involves the square of the nucleon 
and excitation masses. The contribution of excited states for a given tf is smaller for Borel 
improved Sum Rules and a direct comparison of the analysis intervals must take this into 
account. 

One should also acknowledge the presence of an artificial enhancement of the lattice 
correlation functions at very short times. This requires the QCD continuum model to make 
larger contributions relative to the ground state. Any structure in the lattice correlation 
functions at short times reflecting structure in the spectral density is more prevalent in this 
analysis. Hence the enhancement of the correlator gives rise to a more demanding test of 
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the QCD continuum model. Structure not sufficiently accounted for by the model is more 
problematic here. The preceding results give a strong indication that the QCD-SR-inspired 
continuum model provides a sufficiently detailed description of the correlation functions in 
the short time regime to allow the extraction of ground state properties. 

One might consider a similar analysis where the QCD continuum model is replaced by a 
second pole as in (|3.5|) with two poles. Such a comparison would demonstrate the importance 
of the QCD continuum model functional form and the sensitivity of the analysis presented 
here. We consider two fits of (|3.5|) to the lattice data. Since £ was fixed in figures |10| and 



IT| to the optimal value obtained from a four-parameter fit to the time regime 6 — > 20, we 
consider a three-parameter two-pole fit where A^* is similarly fixed. We also present the 
results of a full four-parameter search for comparison. 

For the two-pole spectral density of (|3.5| ), both the ground state nucleon mass and 
coupling strength show a sensitivity to the analysis interval. These parameters do not form 
a plateau when plotted as a function of tf. This failure of the single excitation hypothesis for 
the QCD continuum model is displayed in Figure O for the ground state mass and Figure 
|T5| for the coupling strength. 

Neither the mass nor the coupling strength extracted with a two-pole Ansatz agree with 
the values obtained from the single pole fits to the deep nonperturbative tail of the correlation 
function, where evolution in Euclidean time has isolated the ground state. A successful 
continuum model must reproduce the ground state properties obtained from the simple 
pole fits with tf ~ 20 or 21 ||23|| , as in Figures [10] and [11. For the two-pole spectral density, 



information on ground state properties has been sacrificed in accommodating the deficiencies 
of the spectral density in the short time regime, where the statistical uncertainties are 
relatively small. In this investigation, inclusion of the tail of the correlation function in the 
analysis interval is insufficient to determine the ground state properties when discrepancies 
made by a poor continuum model need to be accommodated. 

In contrast, the QCD-Sum-Rule inspired continuum model produces ground state prop- 
erties that agree with simple pole fits to the tail of the correlation function, and the fit 
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parameters are stable with respect to the analysis interval. Indeed this analysis is very 
sensitive to the continuum model Ansatz. 



VI. COMPARISON WITH OTHER APPROACHES 



A. Operator Product Expansions 



Here we make a comparison of these techniques with those of other approaches for the 
nucleon and A channels. Since QCD Sum Rules and the Random Instanton Liquid Model 
(RILM) [|4[] typically use the interpolator xsr of (|2~3l ) we repeat the previous analysis for 
this interpolator and the A interpolator. The OPEs for these interpolating fields are 

_ , , 3-5/1 2 m„a 1 m 2 ,a 2 



2% 4 V t 6 5 t 5 5 t 4 
47r 2 (: qq :) a 3 



15 t 3 

and 



+ ••• , (6.1) 



^ 3 1 2% 4 U 6 3 t 5 + 3 t 4 



^ (: ^ :)a3 +...), (6.2) 



3 t 3 

where the /x = 2, 3 Lorentz indices of the A interpolator have been summed over. Note that in 
this case the pole contribution to the phenomenological side of (|6.2|) is (4/3) exp(— M&t). 
The QCD continuum models are developed as outlined in Section |T|. A summary of the 
t = 6 — > 20 fit parameters for these interpolating fields is given in Tables y and PTT[ 



B. Systematic Uncertainties 

Prior to comparing with other approaches, it is important to consider possible sources 
of systematic uncertainty, in addition to the statistical uncertainties indicated in the tables. 
As in most exploratory lattice calculations, this investigation considers a single value for the 
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coupling constant (or lattice spacing a) and a single lattice volume. As such, extrapolations 
to the infinite-volume continuum limit are not possible. 

Naively, one might worry that finite volume errors are significant in this investigation. 
The lattice length in the shorter y and z directions is 1.58 fm and is small relative to the 
measured charge diameter, 2 ( r 2 )* , of 1.72(24) fm for the proton |2"5|]2T)|] . However, it is 
important to remember that the actual lattice simulations are performed with quark masses 
larger than the light current quark masses of a few MeV. For the quark masses considered 
on the lattice here, the proton charge diameter ranges from 1.0 to 1.2 fm. For estimating 
finite volume effects, the matter radius, ( (2/3) ( r\ ) + (1/3) ( r\ ) ) 1/,2 ; is a better measure. 
This diameter ranges from 1.0 fm to 1.1 fm, which is about 2/3 the length of the shortest 
dimension of the lattice. 

Order a lattice spacing corrections to the lattice actions can also be a source of systematic 
error. Hadron spectrum calculations indicate these errors are about 10% for hadron 



mass ratios at (3 = 6.0. The main quantity of interest here is the nucleon interpolating field 
coupling strength, X N . Since the nucleon mass is used to determine the lattice spacing, we 
have the advantage of determining the properties of a single hadron. Since M N and 
are linked dynamically, and is constrained to the physical value, we proceed under the 
assumption that finite volume and finite lattice spacing errors in A at are small relative to 
other sources of uncertainty. For the A one should expect systematic uncertainties from 
these effects to be the order of 10%. Quantification of these uncertainties remains a future 
endeavor of lattice QCD investigations. 

The two most prominent sources of systematic uncertainty are the quenched approxima- 
tion and the extrapolation of observables to the chiral limit. Our present understanding of 
the quenched approximation, obtained from numerical simulations of the full theory for a 
few observables, is that it is an excellent approximation to the full theory when the quark 
masses are sufficiently heavy. The quark masses considered in this investigation lie within 
this regime. Hence, it remains to estimate systematic uncertainty in the linear extrapolations 
of masses and coupling strengths. 



Both Mjy and Xn are well behaved in the chiral limit. The leading nonanalytic terms 
of the chiral expansion are nonsingular. Our calculations of the nucleon mass indicate the 
quark mass dependence (or squared pion mass dependence) is linear to a good approxima- 
tion. To order m\, the chiral expansion for is linear in and therefore our lattice 
determinations of the nucleon mass are linearly extrapolated to the chiral limit. Likewise, 
the more phenomenological continuum threshold is taken to have a similar mass dependence. 

To order m 2 , the chiral expansion for A a? is not linear in m 2 . The leading nonanalytic 
term is proportional to m^logm 2 . Here we use a phenomenological description of the 
m 2 behavior of A^ to estimate the size of systematic error in linear extrapolations. For 
the interpolating field xsr it is possible to derive a simple expression for the pion mass 



dependence of A^ [28 



A and c are redundant fit parameters. Here, A is fixed to the inverse lattice spacing. 



Figure [TJ] illustrates the standard linear extrapolation and a fit of Aq and c in fl6.3|) to 



the lattice data. The linear extrapolation and the result from (673) differ by less than one 



standard deviation of the statistical uncertainty in the chiral limit. Hence an estimate of 
the combined statistical and systematic uncertainty for A^v may be obtained by multiplying 
the statistical uncertainties in the tables by 

C. Comparison 

The results summarized in Table |Iy| for Asr evolved to a scale of 1.0 GeV 2 in the leading 
log approximation compare favorably with other approaches utilizing the pole plus QCD 
continuum phenomenology ||29|| . The correct nucleon sum rule predictions of agree 
with the lattice results. Previous works addressing these issues PJ5|]2"^}| refer to the old and 
erroneous results of Ref. || and [[31]]. The numerous errors contained in Ref. are discussed 



and corrected in Ref. 0. The corrections are significant, and restore agreement with these 
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lattice results and with the RILM of \24\. The smaller value of Asr for Ref. || may be due 
to the omission of large a s corrections to the Wilson coefficient of the identity operator used 
to normalize Asr. Such corrections could increase their result by approximately 20%. 

There is some variation in the continuum threshold, s$ from one approach to another. 
This is likely due to differences in the implementation of the QCD continuum models and 
correlation function enhancement in the short time regime reflecting lattice anisotropy. This 
enhancement is expected to induce a larger continuum threshold, and this is confirmed in 
Table [TV|. 

Table [V| also displays uniformity among the residue of the pole for different approaches. 
Here the corrections to Ref. |31| are crucial fl32f . Our analysis displays the usual suppression 



of Ma relative to M^, typical of lattice analyses in which order a and finite volume correc- 
tions remain unaccounted |27|. The large value for Ma from Ref. || is mysterious. Once 



again the value of the continuum threshold is larger than in the QCD Sum Rule analysis as 
anticipated. 

VII. THE PHYSICS OF INTERPOLATORS 

There has been considerable confusion surrounding the physical significance of the baryon 
coupling strength A, and more generally the physics represented in point-to-point correlation 
functions. Here we address these issues. 

In Ref. |24j] some conclusions are drawn from an incorrect interpretation of the physics 



represented in the interpolating field coupling strength, A. From their study of point-to-point 
correlation functions, these authors claim "the octet and decuplet baryons have completely 
different wave functions" and that there is evidence of significant attraction in the scalar 
diquark channel. It must be stated that the phenomenological description of G%(t) in ( |3.12|) 
involves the baryon mass, the continuum threshold describing the effective onset of excited 
states and A, which indicates the ability of the particular interpolator to excite the baryon 
from the QCD vacuum. Information on the nucleon wave function that might be compared 
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with the A wave function is absent in point-to-point correlation functions. 

While it is tempting to compare and A a in hopes of learning something about the 
distribution of quarks in baryons pl |, such a comparison is not correct. It is important to 



recognize that A reflects properties of both the baryon ground state and the interpolating 
field itself. Since the nucleon and A interpolating fields are different it should be no surprise 
to discover A n ^ Aa- Moreover, the nucleon interpolating field is not unique. It is possible 
to construct nucleon interpolators which yield values for A at both larger and smaller than Aa- 
Moreover, their anomalous dimensions are different, making any comparison scale dependent. 
In short, there is little to be learned from a comparison of A a? and Aa- 

Moreover, A does not measure simply the wave function at the origin, but rather the 
wave function at the origin when 

1. the three valence quarks by themselves are in a color singlet state; 

2. the three valence quarks have the particular spin-flavor combination demanded by the 
interpolator. 

The first criterion excludes the majority of the wave function including intermediate states 
of simple diagrams such as a single gluon exchange between quarks. Similarly, the second 
criterion excludes parts of the wave function that have no overlap with the interpolating 
field. These arguments also apply in general to wave function analyses where the locality of 
the baryon annihilation interpolator is relaxed and A is determined as a function of quark 
field operator separations. Of course, wave function analyses also suffer from being either 
dependent on the gauge, or dependent on the path of link variables selected to implement 
gauge invariance. 

The interpolating field acts to bias the physics represented in the wave function. For 
example, if an interpolator with predominant scalar diquark degrees of freedom is used to 
annihilate the baryon, then this part of the full baryon wave function will dominate the 
extracted wave function. The interpolating field filters out the part of the full wave function 
that looks like the interpolator. Thus, wave function analyses allow one to learn about 
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specific sectors of the complete wave function. However, it is dangerous to in turn attribute 
the properties obtained from a specific sector to those of the entire baryon wave function. 

Many have expected it would be possible to find a perfect wave function for the creation 
of a hadron from the vacuum. A perfect smeared source would excite only the ground state 
hadron. Evolution in Euclidean time to suppress excited state contaminations would not 
be required. However, even the best smeared sources require a number of Euclidean time 
steps to isolate the ground state. This inability to obtain a perfect smeared source has a 
simple explanation. The wave function provides information on a particular spin-flavor-color 
Fock-space component of the full wave function. Without information on the full ground 
state wave function, a smeared operator will not be orthogonal to excited states and will 
always generate some excited state contaminations. 

Of course, the only way to probe the properties of hadrons is through the use of a common 
interpolator or current such as the vector current, which has overlap with most hadrons. 
In this way, it is the properties of the hadrons themselves that give rise to differences in 
the extracted observables. For example, vector current matrix elements are determined 
via the calculation of three-point correlation functions in Ref. where a variety 

of electromagnetic observables are reported. These results indicate that scalar diquark 
clustering in the nucleon is actually minimal in QCD p5j|36 . 



VIII. CONCLUSIONS 

A. The QCD Continuum Model 

The physics represented in the short time regime of point-to-point correlation functions 
is described well by the QCD-Sum-Rule-inspired continuum model. The Laplace transform 
of the spectral density appears to be sufficient to render any structure in the spectral den- 
sity insignificant in the short Euclidean time regime of point-to-point correlators. Similar 
conclusions are expected to hold for Borel-improved sum rules. 
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For the lattice correlation functions considered here, both the identity and m q operators 
of the OPE are required when constructing QCD continuum models. Other higher order 
terms of the OPE make negligible contributions to the QCD continuum model due to the 
exponential suppression factor exp(— sot). In general, consideration of the first few leading 
terms of a particular OPE should be sufficient for the construction of the QCD continuum 
model. The QCD continuum model is superior to the use of a second pole. 

B. QCD Sum Rule Test 

The QCD continuum model effectively accounts for excited state contributions to point- 
to-point correlation functions, allowing a determination of ground state properties in the 
short-time regime. However, interplay between pole and continuum model parameters leads 
to rather large uncertainties. 

The smallest fit interval considered here is larger and extends deeper into the nonper- 
turbative regime than that considered in QCD Sum Rule analyses. Hence it would be 
inappropriate to conclude that QCD Sum Rule techniques have been vindicated in this 
analysis. However, the success of the QCD continuum model here is encouraging. This 
analysis justifies the use of the "crude" QCD continuum model and lends credence to the 
results of rigorous QCD Sum Rule analyses. 

C. Future Investigations 

Having established that the QCD continuum model accounts for excited state contribu- 
tions and allows the determination of ground state properties in the short-time regime, future 
lattice calculations may use this technique for analyzing correlation functions that suffer a 
loss of signal prior to a clear ground state domination. In a subsequent paper |J, these 
techniques will be exploited to investigate nucleon properties obtained from unconventional 
interpolating fields. 
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These techniques may be particularly helpful in the analysis of heavy-light meson corre- 
lators. Smeared-source to smeared-sink correlators suffer from large statistical uncertainties, 
while point-to-point correlators suffer from a loss of signal. One of these two correlators are 
needed in order to determine the meson decay constant. In addition, the mass splittings are 
relatively small in these systems, making isolation of the ground state by Euclidean time 
evolution inefficient. On the other hand, smeared-local correlators are well behaved and 
allow a clean determination of the ground state mass. Inclusion of the QCD continuum 
model in the phenomenological description of the correlator allows the analysis of the point- 
to-point correlator and provides a mechanism to account for the problematic excited state 
contaminations. With the pole position previously determined, a reliable extraction of the 
meson decay constant may be possible. 
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TABLES 



TABLE I. ( X1X1 )'• Four-parameter search for the pole plus QCD continuum model. 



Parameter 


n 1 f;o 

K\ — U.lOZ 


, ,. n 1 U 

K2 — U.104 


^3 — U.lOD 


^cr — U.loy o{Z ) 


M N a 

A ia 3 (xl(T 2 ) 

s a 
c 

f from OPE fit 


1.109(8) 

1.17(5) 

1.68(3) 

D.OO^IU J 

5.3(1) 


0.983(8) 

0.94(4) 

1.58(3) 

6.74(9) 

5.6(1) 


0.858(8) 

0.75(3) 

1.49(4) 

6.62(9) 

5.8(1) 


0.628(17) a 
0.38(7) 
1.32(7) 
6.42(19) 


a The physical profc 


on mass sets the lattice 


spacing a = 


0.132(4) fm. 




TABLE II. \ ( 


XSRXsr ) : Four-parame 


ter search for the pole plus QCD continuum model. 


Parameter 


.,. n 1 

K\ — U.lOZ 


k 2 = 0.154 


k 3 = 0.156 


k ct = 0.159 8(2) 


M N a 

Aia 3 (xl0~ 2 ) 
soa 

£ from OPE fit 


1.106(8) 

1.14(5) 

1.63(3) 

5.46(7) 

4.5(5) 


0.980(9) 

0.92(4) 

1.53(3) 

5.25(7) 

4.6(4) 


0.863(9) 

0.76(3) 

1.45(4) 

5.02(7) 

4.6(3) 


0.639(18) 
0.42(8) 
1.29(7) 
4.61(14) 


TABLE III. ( xaXa )'■ Four-parameter search for the pole plus QCD continuum model. 


Parameter 


ki = 0.152 


k 2 = 0.154 


k 3 = 0.156 


k ct = 0.159 8(2) 


M A a 

A A a 3 (xl(T 2 ) 
s a 

£ from OPE fit 


1.164(8) 
2.91(11) 
1.84(2) 
10.3(1) 
7.0L1.5 


1.054(8) 
2.32(9) 
1.78(2) 
10.1(1) 
7.3±1.3 


0.953(9) 
1.82(8) 
1.73(3) 
10.0(1) 
7.6±1.2 


0.758(18) 
0.82(17) 
1.62(5) 
9.8(3) 
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TABLE IV. Comparison with selected results for the nucleoli channel ( XsrXsr )■ 



Approach 


Ref. 


M N 


A SR (l GeV 2 ) 


so 






(GeV) 


(GeV 1 ) 


(GeV) 


This work 




0.96(3) 


0.027(5) 


1.92(11) 


QCD Sum Rules 


Leinweber H 


1.06(18) 


0.031(6) 


1.69(15) 


Lattice (a;-space) 


Chu et al. g 


0.95(5) 


0.022(4) 


<\A 


Instanton Liquid 


Schafer, et al. @ 


0.96(3) 


0.032(1) 


1.92(5) 



TABLE V. Comparison with selected results for the A resonance channel ( xaXa 



Approach 


Reference 


M A 


A A (1 GeV 2 ) 


so 






(GeV) 


(GeV 3 ) 


(GeV) 



This work 
QCD Sum Rules 
Lattice (x-space) 
Instanton Liquid 



Leinweber [32] 
Chu et al. § 
Schafer, et al. (21 



1.13(4) 
1.36 
1.43(8) 
1.44(7) 



0.024(6) 
0.030 
0.037(6) 
0.033(5) 



2.39(10) 
1.58 

3.21(34) 
1.96(10) 
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FIGURES 

FIG. 1. A fit to the lattice correlation function (bullets) from t = 7 — > 20 utilizing a two pole 
Ansatz for the spectral density (solid curve). Individual pole contributions are illustrated by the 
dashed lines. When time slice 6 is included in the fit, the mass of the second pole increases to 
accommodate the earlier time slice. 

FIG. 2. Comparison of the fit of the 1/t 6 term with the ultraviolet cutoff correction (dashed 
line) of ( |3.14 ) and the corresponding curve for 1/t 6 alone (solid line). The lattice cutoff correction 



term is negligible by t — to = 2. In the following, the source and first point of the correlator are 
discarded and all fits begin at t = 6. 

FIG. 3. QCD continuum model contributions. From top down the lines represent the total 
continuum model contribution, followed by the individual contributions of the identity operator, 
m q correction, quark condensate and finally the ml correction. 

FIG. 4. OPE contributions. From top down, at time slice 6, the lines represent the sum 
total of the first four terms of the OPE, followed by the individual contributions of the identity 
operator, m q correction, quark condensate and finally the mi term. Line types are as in Figure ||. 
The crossing of the quark condensate contribution with the identity operator at t = 8 indicates the 
breakdown of the truncated OPE. The shaded areas indicate the regions typically excluded from 
Borel-transformed QCD Sum Rule analyses. 



FIG. 5. The two-point correlator at k = 0.152 for the nucleon interpolating field xi of fl2,l|) . 
The fits for the 12 analysis intervals are illustrated. The source position is at to = 4. Neither the 
source nor t = 5 are included in the fit as indicated in the discussion surrounding Figure |2[ 

FIG. 6. The nucleon mass determined in each analysis interval plotted as a function of tf. In 
this and the following figure, bullets correspond to pole plus QCD continuum fits from t = 6 — > tf, 
and open squares, offset for clarity, illustrate a simple pole fit to the region tf — 7 — > tf which was 
selected to give similar statistical uncertainties in the nucleon mass at t = 20. 
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FIG. 7. The nucleon coupling strength determined in each analysis interval plotted as a 
function of tf. Symbols are as in Figure ||. 

FIG. 8. The quark mass dependence of £ for mean- field-improved wave function renormaliza- 
tion (solid curve) and naive wave function normalization (dashed curve). 

FIG. 9. Correlation function fits for each of the 15 analysis intervals t = 6 — > tf, where 
t f = 9 -> 23 for the QCD-SR method test. 

FIG. 10. The nucleon mass determined in each analysis interval plotted as a function of tf. In 
this and the following figure, bullets correspond to pole plus QCD continuum fits from t = 6 — ► tf, 
and open squares illustrate a simple pole fit to the region tf — 7 — > tf as in Figure |6|. 

FIG. 11. The nucleon coupling strength determined in each analysis interval plotted as a 
function of tf. Symbols are as in Figure [0]. 

FIG. 12. The nucleon ground state mass determined in each analysis interval, t = 6 — ► tf, 
plotted as a function of tf for a spectral density involving two poles, as opposed to a pole plus 
QCD continuum model. Bullets correspond to a three-parameter fit with Aat* fixed as described 
in the text. Open circles illustrate the full four-parameter fit. The extracted mass depends on the 
analysis interval and fails to agree with the simple pole fits (open squares) from tf — 7 — » at 
tf ~ 20 or 21. This confirms the sensitivity of this analysis to the QCD continuum model Ansatz. 

FIG. 13. The nucleon coupling strength plotted as a function of tf for a spectral density 
involving two poles. Symbols are as in Fig. [12|. The dependence on the analysis interval and failure 
to agree with the simple pole fits (open squares) at tf ~ 20 or 21 confirms the sensitivity of this 
analysis. 

FIG. 14. Extrapolation of Asr to the chiral limit. The solid line indicates the standard linear 
extrapolation, while the dashed curve illustrates the two-parameter fit of fl6.3| ) to the lattice data. 
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